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Abstract Magnetoelectronics is mainly digital, i.e. governed by up and 
down magnetizations, hi contrast, analogue magnetoelectronics makes use 
of phenomena occuring for non-coUinear magnetization configurations. Here 
we review theories which have recently been applied to the transport in non- 
collinear magnetic nanostructures in two and multiterminal structures, viz. 
random matrix and circuit theory. Both are not valid for highly transparent 
systems in a resistive environment like perpendicular metallic spin valves. 
The solution to this problem is a renormalization of the conventional and 
spin-mixing conductance parameters. 



1 Introduction 

The giant magnetoresistance, as well as most of the current magnetoelec- 
tronics, can be understood in terms of the transport of electrons in either 
spin-up or spin-down state, since the magnetizations are collinear (paral- 
lel or antiparallel) with the spin-quantization axis. Both charge and spin 
transport can be described in terms of two "channels" with spin-dependent 
conductivities, scattering rates etc. fD. This "digital" magnetoelectronics 
does not profit from the "analogue" freedom of a magnetization to point in 
any direction. Early seminal contributions by Slonczewski |21 and Berger P] 
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revealed fundamentally new physics and technological possibilities of non- 
collinearity, which triggered a large number of experimental and theoretical 
studies. An important example is the non-equilibrium spin-current induced 
torque (briefly, spin torque) which one fcrromagnet can exert on the mag- 
netization vector of a second magnet through a normal metal in a biased 
spin valve structure. This torque can be large enough to dynamically turn 
magnetizations ^ , which is potentially interesting as a low-power switching 
mechanism for magnetic random access memories 0. Non-collinear magne- 
tizations are also essential for novel magnetic devices like the spin-flip jSlE] 
and spin-torque transistors detection of spin-precession the Gilbert 
damping of the magnetization dynamics in thin magnetic films |10|. and 
spin-injection induced by ferromagnetic resonance 

We are interested in heterostructures containing band ferromagnets that 
are accurately described by a Stoner spin-density functional model. Elemen- 
tal metals and its alloys have high electron densities and their thin-film het- 
erostructures are usually considerably disordered. Size quantization effects 
on transport can therefore mostly be disregarded Semiclassical meth- 
ods are appropriate in slowly varying bulk regions of the structures, but 
heterointerfaces can be atomically sharp and must be treated fully quan- 
tum mechanically. There are basically two methods which are suitable to 
understand and compute transport properties of these systems from first 
principles, viz. Green function theory with configurational averaging 
and the scattering formalism for transport, combined with random matrix 
theory These two approaches have recently been extended to non- 
collinear magnetic structures, viz. magnetoelectronic circuit ;6 and random 
matrix theory '15' . Here, as in '16', we show that both approaches are closely 
related, but do not hold for transparent interfaces. Following Schep's [T7| 
strategy for coUinear systems, both theories can be generalized, leading to 
analytical results for perpendicular spin valves with parameters that can 
be obtained by ab initio band structure calculations, as well as determined 
quantitatively by experiments. 

2 Boltzmann and Diffusion Equation 

When a local non-equilibrium magnetization docs not point in the direction 
of the spin-quantization axis, the distribution function for band states at 
the Fermi energy with index n is a matrix in Pauli spin space 

On the right hand side the distribution is expanded into unit matrix and 
the vector of the Pauli spin matrices. is charge accumulation and the 
spin accumulation ff^ is a vector whose direction is always parallel to the 
magnetization vector m in the bulk of a ferromagnet, but arbitrary in a 
normal metal depending on device configuration and applied biases. /„ can 



Digital to analogue 



3 



be diagonalized by unitary rotation matrices in spin space, characterized 
by the polar angles 6 and ip. Let us assume for simplicity that these angles 
are piecewise constant in position space, thus disregarding magnetic domain 
walls proi and magnetic field- induced spin precession in normal metals 
In the local spin quantization frame the distribution function is diagonal 
with two spin components s = ±1. Introducing spin-conserving and spin- 
flip scattering life times Ts and t|^_3, for the sake of argument taken to be 
state- independent, and separating the distribution into an isotropic elec- 
trochemical potential /i^ and an anisotropic term 7„s that vanishes when 
averaged over the Fermi surface, the Boltzmann equation for the stationary 
state reads [T^ 

Ts,-s J '''s,-s 

where v„s is the group velocity of state ns. Charge and spin currents read 

jc^ -^^^nslns , (3) 

ns 

ns 

The Boltzmann equation is still unnecessarily complicated for most realistic 
systems. In the presence of sufficient disorder, only the lowest harmonics 
of 7„s in reciprocal space survive. In that limit, the Boltzmann equation 
reduces to the diffusion equation 



sd 



= y/Dr^ is the spin-flip diffusion length, which does not depend on spin 
index I2Q| . The spin-averaged diffusion coefficient D can be written in terms 
of the density of states at the Fermi energy Ns [Ep) 

{N^{Ef)+Ni{Ef))D N^{Ef)D^ Ni{Ef)Di' 

in terms of the spin-dependent diffusion coefficients. In a simple two-band 
model Ds = VsTs/<i, where Vs are the spin-dependent Fermi velocities. The 
average spin-flip relaxation time is defined as 

The currents 

Ur) = -^VM,(r) (8) 
are governed by the spin-dependent conductivities 

<Js^ NsiEF)e^Ds . (9) 
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Fig. 1 Different realizations of perpendicular spin valves in which 6 is the angle 
between magnetization directions, (a) Highly resistive junctions like point contacts 
and tunneling barriers limit the conductance, (b) Spin valve in a geometrical con- 
striction amenable to the scattering theory of transport, (c) Magnetic multilayers 
with transparent interfaces 



3 Boundary Conditions 

The semiclassical approach is valid when the potential landscape varies 
slowly on the scale of the Fermi wave length. In heterostructures we of- 
ten encounter regions in which materials change on an atomic scale, such 
as at intermetallic interfaces or tunnel junctions, which have to be treated 
quantum mechanically. The "nodes" are the bulk regions, in which the semi- 
classical distributions are well defined. The intermediate scattering regions, 
or "contacts", can then be treated formally exactly by boundary condi- 
tions, which link the distributions of two neighboring nodes. Consider the 
spin valve structures in Fig. 1, which may be part of a larger circuit. We de- 
note the distribution functions in the ferromagnetic terminals by subscripts 
L end R. We explicitly allow for a drifting distribution by the superscript 
a = ±1, which indicates whether drift is in (a = 1) or against (a = — 1) 
the transport direction (from left to right). Taking into account the differ- 
ence between the left- and right-moving distribution functions is the key 
generalization of the previous theories in our treatment here. In Sec. 5 it 
is discussed how the circuit theory can be recovered be renormalizing the 
conductance parameters. Here we concentrate on scattering theory. 

In order to work with simple matrices instead of diadics, we follow Wain- 
tal et al. |15j and introduce the 4x1 vector representation [f" (r)]"^ — 
(/TTW)/n(i")'/Ti(r)>/ii (!■))"• The boundary conditions for the non- 
equilibrium distributions to the left and right of the scattering region then 
read: 

f^n = Yl (^L^") n,n ^t™ + (^^^r) nm ^R.m ' (10) 

meL meR 

fZ.n = Y (^i^i)nm ^L,m + Y i^"'-*^) .^„^ ffl.,„ • (H) 
meL jeR 

T, R are 4x4 transmission and reflection probability matrices and the sub- 
scripts indicate the direction of the currents {L —t R denotes transmission 
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from left to right, i? — > i? reflection from the right, etc.). All matrix ele- 
ments follow from the scattering matrix and are conveniently normalized, 
for example 

The transmission amplitudes, such as t^^^s' of a wave coming in from 
the left as mode m and spin s' and going out in mode n and spin s, are 

here collected in vectors t^;^^ (t^f'^',tff"',tf^^,tff^'j . We also 

have S e [1, 4], iVl = Nf {Ss.i + 63.2) + N[ {5s,3 + 63,4) , where Vf is the 
number of modes for spin s in the ferromagnet. 

In a nutshell, this is a very general formulation of charge and spin trans- 
port, but it is not yet amenable for analytic treatment or analysis of exper- 
iments. The isotropy assumption that reduced the Boltzmann to the diffu- 
sion equation in the previous chapter, enormously simplifies the results, as 
demonstrated in the following. 

We focus here on the electrical charge current as a function of the magne- 
tization configuration in symmetric spin valves, as in Fig. 1(b), (c), in order 
to keep the analytical manipulations manageable. We will see later that we 
can derive rules from these results that are valid for general structures. T 
and R are functions of the magnetic configuration, which, disregarding mag- 
netic anisotropics, can be parameterized by a single polar angle 9. In first 
instance, we disregard spin-flip scattering and discuss later how it can be in- 
cluded. Integrating over the lateral coordinates leaves a position dependence 
only in the transport direction (x). The next step is the assumption that 
the distribution functions for incident electrons from the left and right are 
isotropic in space. The distribution functions for the outgoing electrons do 
not have to be isotropic, as long as they are subsequently scrambled in the 
nodes. The isotropy assumption may be invoked when the nodes are diffuse 
or chaotic, such that electrons are distributed equally over all states at the 
(spin-dependent) Fermi surfaces (which is equivalent to replacing state de- 
pendent scattering matrix elements by its average ^3). The Fermi surface 
integration is then carried out easily, and the distribution functions within 
left and right ferromagnet nodes (at locations xl and xr, respectively) are 
matched via simplified boundary conditions 

f+ {xr) = Tl^r {9) f+ [xl] + Rr^r (9) i- (xr) , (13a) 
f " (xl) = Rl^l (9) f+ (xl) + Tr^l (9) i- (xr) , (13b) 

where the 4x4 transmission and reflection probability matrices have ele- 
ments like ^ni^ 

S mn 
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In the coordinate systems defined by the magnetization directions, the trans- 
verse components of the spin accumulation in the ferromagnets vanish iden- 
tically [Hll21| and the distributions in the magnets depend on the local spin 
current densities 7s and chemical potentials /is only 

K/i?) = ((±7T + /^t) K/fl) ' 0' 0' (±7i + Mi) K/fl)) . (15) 

By this choice the explicit angle-dependence of transport is contained only 
in the matrices. 

We can now link an arbitrary distribution on the left to compute the 
distributions on the right, subject to the constraint of charge current con- 
servation. Here we focus on the simple case in which we apply a bias 

s 

but no spin accumulation gradient /is (xl) — ^J--s (xl) — fJ-s (xr) — ^i-s {xr) 
over the system. We then find that 7s (2;^) = 7s {xr) , i.e. the spin current 
component parallel to the magnetization on left and right ferromagnets are 
the same. The charge current 

divided by the chemical potential drop is the electrical conductance G = 
Ic/A^. Eqs. 113I15|) then lead to 

G = ¥ E [i - Tl-,r + Rr^r] Tl^r} . (18) 

5=1,4 -^-^ 
S'=1A 

When the transparency is small, all transmission probabilities are close to 
zero, reflection probabilities are close to unity, and the Landauer-Biittiker 
conductance, starting point of is recovered: 

S=l,4 
S' = l,4 

Indeed, in this limit the distributions to the left and right are not perturbed 
by the current, the nodes are genuine reservoirs, and standard scattering 
theory applies. Also, when 9 = 0,7r, Eq. l{TH|l is equivalent to results by 
Schep et al. ^T] for the two-channel model. 

The scattering region is still not specified and may be interacting and/or 
quantum coherent. We now discuss how analytical results can be obtained 
in the non- interacting, diffuse limit. 
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4 Semiclassical Concatenation 

The scattering matrix of a composite system can be formulated as concate- 
nations of the scattering matrix from separate elements, e.g.. the scatter- 
ing matrices of bulk layers and interfaces |22| . By assuming isotropy, i.e. 
sufficient disorder or chaotic scattering, Waintal et al. ^Sj proved by av- 
eraging over random scattering matrices that size quantization effects like 
the equiUbrium exchange coupling or other phase coherent phenomena are 
destroyed by disorder and vanish like the inverse of the number of modes. 
Under these conditions we are free to define nodes in the interior of the de- 
vice and link them via the boundary conditions H13|) . This is equivalent to 
composing the total transport probability matrices in Eq. H18() in terms of 
those of individual elements by semiclassical concatenation rules jSJ- The 
4x4 transmission probability matrix through a F(0)/N/F(6') double het- 
erojunction as in Fig. 1 in which bulk scattering is absent, takes the form 

f [9) = Tn^f {9) [I - i?N^N (0) i?N^N {0)] Tp^N (0) , (20) 
where the interface transmission and reflection matrices as a function of 
magnetization angle appear. The transformations needed to obtain Tn^f [9) 
and i?N^N (^) require some attention. In terms of the spin-rotation 

~ /cos 0/2 -sin 0/2 \ 
^ \ sui9/2 cos6'/2 ) ^ ' 

and projection matrices (s = ±1) 

u.(^)=-f'+'T^ ''"\) ' (22) 

^ ' 2 \ s sm 1 - s COS y ' ^ ' 

the interface scattering coefficients (omitting the mode indices for simplic- 
ity) are transformed as follows 6 

Ut^mU^ ^Y.^^^'s'' , (23) 



rN^N 



ti:r^ = UsAf , (24) 

t^r^ = t^Ul , (25) 

rf,r^ = rfSss' . (26) 

The superscript c indicates that the matrices should be evaluated in the ref- 
erence frame of the local magnetization, and are thus diagonal in the absence 
of spin-flip relaxation scattering at the interfaces. Different transformation 
properties for the different elements of the scattering matrix derive from our 
choice to use local spin-coordinate systems that may differ for each magnet. 
Let us, for example, inspect a transmission matrix element from the normal 
metal into the ferromagnet with magnetization rotated by 9 



T mn T 



(27) 
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and analogously for the other matrix elements as well as other matrices. 

Transport through a more complex system can be treated by repeated 
concatenation of two scattering elements in terms of reflection and trans- 
mission matrices analogous to Eq. H2U|I . In the presence of significant bulk 
scattering, we can represent a disordered metal B with thickness cIb by 
diagonal matrices like ;17iil5) 

where pf , Ab are the single-spin bulk resistivities and cross section of the 
bulk metal (normal or magnetic). 

The interface parameters of the present theory are the spin-dependent 
Landauer-Biittiker conductances 



Im.s 



hn Irn hn Im 



(29) 

and the real and imaginary part of the spin-mixing conductance 

Im 

which can also be represented in terms of the total conductance <? = 5t + ; 
polarization p ^ [g^ — g{) / g, and relative mixing conductance rj = ^g^^j g. 

The actual concatenation of the 4x4 matrices defined here is rather 
complicated even when using symbolic programming routines. This explains 
why in |15| analytic results were obtained only in special limiting cases. We 
found that final results are simple even in the most general cases, not only 
for Eq. (|19|l considered by ^S], but also for Eq. (|18|l . For the spin valves in 
Fig. 1, we find for the conductance as a function of angle 



where f] = 2g^y/g and 



1 


1 


PF,sdF 


9s 




h 2Af 


1 


1 


PNdN 


9U 


9n 


h 2An 



Re fi 1 — cos 6 



PNdN 1/1,1 



h 2Af 2 \Nf Nn 
1 



2Nn 



(31) 
(32) 



Equation H3U|) is identical to the angular magnetoresistance derived by cir- 
cuit theory 6 after replacement of (js and by gs and g-^^. Physically, 
in Eqs. (|31I32|I spurious Sharvin resistances are substracted from the inter- 
face resistances obtained by scattering theory, whereas bulk resistances are 
added. These corrections are large for transparent interfaces and essential 
to obtain agreement between experimental results of transport experiments 
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in CPP (current perpendicular to plane) multilayers |25II26II27| and first- 
principles calculations, for conventional J17','28"'29" as well as mixing con- 
ductances .7,. The mixing conductance parameterizes the magnetization 
torque due to a spin accumulation in the normal metal, governed by the re- 
flection of electrons from the normal metal. It is therefore natural that the 
mixing conductance is reduced by the bulk resistance of the normal metal 
and we can also understand that only the normal metal Sharvin resistance 
has to be substracted. The real part of the mixing conductances is often 
close to the number of modes in the normal metal gi^i « A^at, in which case 
g^l w Pni- By letting Nf ^ oo we are in the regime of ^B]. The 

circuit theory is recovered when, additionally, Nn — > oo. The bare mixing 
conductance is bounded not only from below Reg^^ ^5/2 [EIj but also from 
above \guf /^'^9n ^ 2iVjv 



5 Extended Circuit Theory 

It is not obvious how these results should be generalized to more compli- 
cated circuits and devices as well as to the presence of spin-flip scattering in 
the normal metal. The magnetoelectronic circuit theory [Sj does not suffer 
from these drawbacks. Originally, it was assumed in Ref. |H] that local spin 
and charge currents through the contacts only depend on the generalized 
potential differences, and the local node chemical potentials are obtained 
by a spin-generalization of the Kirchhoff laws of electrical circuits. This is 
valid only for highly resistive contacts, such that the in and outgoing cur- 
rents do not significantly disturb the quasi-equilibrium distribution of the 
nodes. Fortunately we are able to relax this limitation and take into account 
a drift term in the nodes as well. In order to demonstrate this, we construct 
the fictitious circuit depicted in Fig. 2. Consider a junction that in conven- 
tional circuit theory is characterized by a matrix conductance g, leading 
to a matrix current i when the normal and ferromagnetic distributions /l 
and fii are not equal. When the distributions of the nodes are isotropic, we 
know from circuit theory that 

i = ^9)ss' Us (/l - /fl) Us' , (33) 

ss' 

where the projection matrices Us are defined in Eq. 122|1 and {g)ss — 9s, 
{g)s.-s =9s.-s- Introducing lead conductances, which modify the distribu- 
tions /z, — > fi and /2 <— Jr, respectively, we may define a (renormalized) 
conductance matrix g, which causes an identical current i for the reduced 
(matrix) potential drop: 

*' = ^(g)^^,iis(/i-/2)^is'. (34) 

ss' 
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Fig. 2 Fictitious device that illustrates the generalization of circuit theory to 
transparent resistive elements as discussed in the text 



When the lead conductances are now chosen to be twice the Sharvin con- 
ductances, and using (matrix) current conservation 



(35) 
(36) 



= ^2<u, (/2-/fl,) 



straightforward matrix algebra leads to the result that the elements of 
g are identical to the renormalized interface conductances found above 
[Eqs. I|31I32|I without the bulk resistivities]. By replacing g by g we not 
only recover results for the spin valve obtained above, but we can now use 
the renormalized parameters also for circuits with arbitrary complexity and 
transparency of the contacts. Also spin-flip scattering in N can now also be 
included 0; it does not affect the form of Eq. (|37|l either, but only reduces 
the parameter Other effects of the spin-flip scattering arc discussed in 
detail by Kovalev et al. 



6 Applications 

Intermetallic interfaces in a diffuse environment (see Fig. Ic) have been stud- 
ied thoroughly by the Michigan State University collaboration jJS] and oth- 
ers |26II27| in perpendicular (CPP) spin valves. These experiments provided 
a large body of evidence for the two-channel {i.e. spin-up and spin-down) 
series resistor model and a wealth of accurate transport parameters such as 
the spin-dependent interface resistances for various material combinations 
|17II27II28II29| . In exchange-biased spin valves, it is possible to measure the 
electric resistance as a function of the angle between magnetizations, which 
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has been analyzed experimentally and theoretically |81IIH2j . Pratt c.s. ob- 
served that experimental magnetoresistance curves \^^'^\ could accurately be 
fitted by the form ^ 

i?(7r)-i?(0) x(l + cose')+2 ^ ' 

According to the new insights described above, the free parameter x is a 
function of renormalized microscopic parameters 

X=r^]fC-l (38) 

in terms of the relative mixing conductance rj — 2giji/g, the polarization 
P = (.91 ~ 9 1 ) /S' ^^'^ average conductance g = + gi- 

Experimental values for the parameters for Cu/Permalloy (Py) spin 
valves are x — 1-2 and p — 0.6 . Disregarding a very small imagi- 
nary component of the mixing conductance 0, using the known values 
for the bulk resistivities, the theoretical Sharvin conductance for Cu (0.55 • 
10^^ ri^^m^^/spin 17 ), and the spin-flip length of Py as the effective thick- 
ness of the ferromagnet {£f^ — 5 nm j25j), we arrive at the bare Cu/Py in- 
terface mixing conductance Gn = 0.39 (3) •10^'^ il~^m~^. This value may be 
compared with the calculated mixing conductance for a disordered Co/Cu 
interface (0.55 • IQ-*^^ r2~-'^m [7]). The agreement is reasonable, but leaves 
some room for material and device dependence that deserves to be investi- 
gated in the future. The mixing conductance can also be determined from 
the excess broadening of ferromagnetic resonance spectra. A larger mixing 
conductance in Pt /Py can be explained by the larger density of conduction 
electrons in Pt compared to Cu (201 . Reasonable agreement between exper- 
iment and theory has been also found by Zwierzycki et al. |34j for Fe/Au 
multilayers. 

The spin torque on a ferromagnet pill5| equals the spin current through 
the interface with vector component normal to the magnetization direction 
and its evaluation is closely related to the charge conductance |51ll5| . An 
analytical expression for the spin valve reads: 

L{e) = ^-^^^ ^ (39) 

l-cos6i+J^(l + cos6') 87r 

Note that here the imaginary part of the mixing conductance is taken into 
account explicitly, but the torque remains coplanar to the magnetization 
of the contacts, i.e. an out-of-plane "effective" field vanishes identically. 
Previous results pill5| are recovered in the limit that fj ^ 2 and p ^ 1. 
By the generalized circuit theory it is now straightforward to compute the 
torque on the base contact of the spin-flip transistor with antiparallel source- 
drain magnetizations [7] . Let us assume that the three contacts are identical, 
and the base contact magnetization lies in the plane of the source and drain 
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Fig. 3 The spin-accumulation induced magnetization torque for a two-terminal 
spin valve and a three-terminal spin-flip transistor. Ajj, is the source-drain bias 
and all contact parameters are taken to be the same, with Re rj — 2 and Im = 



magnetizations. Assuming that we may disregard spin-flip scattering in the 
base contact, the in-plane torque Lb turns out to be always larger than 
the spin valve torque L in the two-terminal spin valve, Eq. H39(l . with a 
symmetric and flatter dependence on the angle of the base magnetization 
direction 6 (see Fig. 3) 

L (0)= gpRefisind Afi 

(l-Re?7)cos2 6l-f Rer)-|-— -^f— — T 47r' 

^ ' 2+|r;| /(R,e»;) 

In the presence of a significant imaginary part of the mixing conductance, 
we also find an out-of-plane (effective field) torque L± [0) with the same 
angular dependence and 

Lb \r^f + 2Yiefi' 

Stiles and Zangwill ^35. directly solved the Boltzmann equation for spin 
valves to obtain angular magnetoresistance and spin torque, approximat- 
ing the mixing conductance by the number of modes (note that in a direct 
solution of the Boltzmann equation this parameters should not be renormal- 
ized). The numerical results agree well with the functional form (|37|l (M.D. 
Stiles, private communication) . This function has been also derived by Slon- 
czewski j^E) and later by Shpiro et al. |SZ| ■ Sloncewski rederived this result 
with a simple circuit theory similar to that of and also pointed out the 
relation between the angular magnetoresistance and the spin torque. Shpiro 
et al. |37] found the form ()37|l to be valid in the limit of vanishing exchange 
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splitting, thus in a regime different from the transition metal ferromagnets 
considered here PT] , 

7 Conclusions 

We reported analytical results for the angular magnetoresistance of arbi- 
trary spin valves, which, by comparison with experiments |33| . leads to a 
value for the mixing conductance and spin torque for the Cu/Py interface 
of G-^i — 0.39(3) • 10^^ il^^m^^. The associated generalization of magne- 
toelectronic circuit theory opens the way to engineer materials and device 
configurations to optimize switching properties of magnetic random access 
memories. Mixing conductances determined by experiments or first princi- 
ples theory are transferable to arbitrary devices and may be used for static 
as well as dynamic transport properties. The spin-dependent interface re- 
sistances determined by CPP-GMR transport experiments have played an 
important role in understanding "digital magnetoelectronics" . We hope that 
the spin-mixing conductances will play a comparable role in "analogue mag- 
netoelectronics" . 

Acknovifledgements We profited from discussions with Bart van Wees, 
Paul Kelly, Alex Kovalev, and Yuli Nazarov, and have been supported by 
FOM, NSF Grant DMR 02-33773 and the NEDO joint research program 
"Nano-Scale Magnetoelectronics" . 

References 

1. S. Maekawa and T. Shinjo (eds.), Spin Dependent Transport in Magnetic 
Nano structures (Taylor and Francis, London, 2002). 

2. J.C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996). 

3. L. Berger, Phys. Rev. B 54, 9353 (1996). 

4. M. Tsoi, A. G. M. Jansen, J. Bass, W.-C. Chiang, M. Seek, V. Tsoi, and 
P. Wyder, Phys. Rev. Lett. 80, 4281 (1998); J.-E. Wegrowe, D. Kelly , Y. 
Jaccard , P. Guittienne, and J.-P. Ansermet, Europhys. Lett. 45, 626 (1999); 
J.Z. Sun, J. Magn. Magn. Mater. 202, 157 (1999); E.B. Myers, D.C. Ralph 
, J.A. Katine, R.N. Louie, and R.A. Buhrman, Science 285, 867(1999); J.A. 
Katine, F.J. Albert, R.A. Buhrman, E.B. Myers, and D.C. Ralph , Phys. 
Rev. Lett. 84, 3149 (2000); J. GroUier, V. Cros, A. Hamzic, J.M. George, H. 
Jaflres, A. Pert, G. Faini, J. Ben Youssef, and H. Legall, Appl. Phys. Lett. 
78, 3663 (2001); B. Oezyilmaz, A.D. Kent, D. Monsma, J.Z. Sun, M.J. Rooks, 
and R.H. Koch, cond-mat/0301324 S. Urazhdin, N.O. Birge, W.P. Pratt Jr, 
and J. Bass, cond-mat/0303149 . 

5. K. Inomata, lEICE Transactions on Electronics E84-C, 740 (2001); J.C. Slon- 
czewski, cond-mat/0205055 

6. A. Brataas, Yu.V. Nazarov, and G.E.W. Bauer, Phys. Rev. Lett. 84, 2481 
(2000); Eur. Phys. J. B 22, 99 (2001). 

7. K. Xia, P. J. Kelly, G.E.W. Bauer, A. Brataas, and L Turek, Phys. Rev. B 
65, 220401 (2002). 



14 



Gerrit Bauer et al. 



8. G.E.W. Bauer, A. Brataas, Y. Tserkovnyak, and B.L van Wees, Appl. Phys. 
Lett., in press. 

9. D. Huertas-Hernando, Yu.V. Nazarov, A. Brataas, and G.E.W. Bauer, Phys. 
Rev. B 62, 5700 (2000). 

10. Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, Phys. Rev. Lett. 88, 117601 
(2002). 

11. A. Brataas, Y. Tserkovnyak, G.E.W. Bauer, and B. Halperin, Phys. Rev. B 
66, 060404 (2002). 

12. The quantum size effects observed by S. Yuasa, T. Nagahama, and Y. Suzuki, 
Science 297, 234 (2002) are an exception, because the high-quahty tunnel 
barriers focus conductance electrons to essentially a single wave vector, so the 
Fermi surface averaging that smears out residual quantum effects in metallic 
systems, is not effective. 

13. J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

14. C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

15. X. Waintal, E.B. Myers, P.W. Brouwer, and D.C. Ralph, Phys. Rev. B 62, 
12317 (2000). 

16. G.E.W. Bauer, Y. Tserkovnyak, D. Huertas-Hernando, and A. Brataas, Phys. 
Rev. B 67, 094421 (2003), cond- mat/0205453 

17. K.M. Schep, J.B.A.N. vanHoof, P.J. Kelly, G.E.W. Bauer, and J.E. Inglesfield, 
Phys. Rev. B 56, 10805(1997); G.E.W. Bauer, K.M. Schep, P.J. Kelly, and 
K. Xia, J. Phys. D: Appl. Phys. 35, 2410 (2002). 

18. G. Tatara and H. Fukuyama, Phys. Rev. Lett. 78, 3773 (1997). 

19. T. Valet and A. Pert, Phys. Rev. B 48, 7099 (1993). 

20. A. Filip, Ph.D. Thesis, University of Groningen, 2002. 

21. M.D. Stiles and A. Zangwill, Phys. Rev. B 66, 014407 (2002). 

22. S. Datta, Quantum Phenomena in Semiconductor Nanostructures, (Addison 
Wesley, 1989). 

23. B. Shapiro, Phys. Rev. B 35, 8256 (1987); M. Cahay, M. McLennan, and S. 
Datta, Phys. Rev. B 37, 10125 (1988); A. Brataas and G.E.W. Bauer, Phys. 
Rev. B 49, 14684(1994). 

24. A.A. Kovalev, A. Brataas, and G.E.W. Bauer, Phys. Rev. B 66, 224424 
(2002). 

25. W.P. Pratt, Jr., S.-F. Lee, J.M. Slaughter, R. Loloee, P.A. Schroeder, and J. 
Bass, Phys. Rev. Lett. 66, 3060 (1991); J. Bass and W.P. Pratt, J. Magn. 
Magn. Mater. 200, 274 (1999). 

26. M.A.M. Gijs, S.K.J. Lenczowski, and J.B. Giesbers, Phys. Rev. Lett. 70, 3343 
(1993). 

27. For recent reviews see: M.A.M. Gijs and G.E.W. Bauer, Advances in Physics 
46, 285 (1997); J.-P. Ansermet, J. Phys.-Cond Mat. 10, 6027 (1998); A. 
Barthelemy, A. Fert and F. Petroff, in Handbook of Magnetic Materials, Vol. 
12, edited by K.H.J. Buschow (1999); E. Tsymbal and D. G. Pettifor, Sol. 
State Phys. 56, 113 (2001). 

28. M. D. Stiles and D. R. Penn, Phys. Rev. B 61, 3200 (2000). 

29. K. Xia, P.J. Kelly, G.E.W. Bauer, I. Turek, J. Kudrnovsky, and V. Drchal, 
Phys. Rev. B 63, 064407 (2001). 

30. Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, Phys. Rev. B 66, 224403 
(2002). 

31. P. Dauguet, P. Gandit, J. Chaussy, S.F. Lee, A. Fert, and P. Holody, Phys. 
Rev. B 54, 1083 (1996). 



Digital to analogue 



15 



32. A. Vedyayev, N. Ryzhanova, B. Dieny, P. Dauguet, P. Gandit, and J. Chaussy, 
Phys. Rev. B 55, 3728 (1997). 

33. L. Giacomoni, B. Dieny, W.P. Pratt, Jr., R. Loloee and M. Tsoi, to be pub- 
lished. 

34. M. Zwierzycki et al., unpublished. 

35. M. D. Stiles and D. R. Penn, J. Appl. Phys. 91, 6812 (2002). 

36. J.C. Slonczewski, J. Magn. Magn. Mat. 247, 324 (2002). 

37. A. Shpiro, P.M. Levy, and S. Zhang, cond-mat /0212045. 



